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Abstract. We prove the global-in-time existence of weak solutions to the Navier-Stokes equations 
of compressible isentropic flow in three space dimensions with adiabatic exponent 7 > 1. Initial 
data and solutions are small in around a non-constant steady state with densities being positive 
and essentially bounded. No smallness assumption is imposed on the external forces when 7 = 1. 
A great deal of information about partial regularity and large-time behavior is obtained. 
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1. Introduction 

We prove the global existence of weak solutions to the Navier-Stokes equations of compressible flow 
in three space dimensions: 

f pt+div{pu) =0, , 
\ {pu^)t + dW{pu^u) + (P^),^ - /i Au -I- (e - p) (div u),^ +p. ^ ■ ' 

Here p and u = {u^,u'^, u^) are the unknown functions of a; G M'^ and t > 0, P-^ — P-y{p) is the pressure, 
/ is the external force, p and S, — p are viscosity constants. 
The system (|l.ip is solved subject to initial conditions 

(p(-,0),u(-,0)) = (po,wo), (1.2) 

where po is bounded above and below away from zero, uq G for some p > 6, and modulo constants, 
(po ~ PijUi)) is small in L'^iV?) for some specific non-constant function p~^ which will be defined 
later. The solvability to various Cauchy problems for the Naiver-Stokes equations has been discussed 
by many other mathematicians for decades. Matsumura-Nishida 9 proved the global existence of 
solutions around a constant state when the initial data is taken to be small in , and later 
Danchin [1 generalized their results by replacing with certain Besov spaces of functions. On the 
other hand. Lions 8 and Feireisl [2]-[3] proved global existence of weak solutions to (|l.ip - (|1.2|) with 
finite-energy initial data and nonnegative initial density. In between those two type of solutions as 
mentioned, Hoff [4j-[6j developed a new theory of intermediate regularity class solutions which may 
exhibit discontinuities in density and velocity gradient across hypersurfaces in R'^. In the presence of 
large external force, our work generalizes earlier results of Matsumura-Yamagata [IQ^ in two ways: the 
restriction on the L°° norm oi p — p^ has been eliminated and initial velocity is not necessary in . 

We introduce two variables associated with the system (jl.ip which are important to our analysis. 
The first one is the usual vorticity matrix uj ~ u^'^ — u^^ — u*^. , while the other one is the effective 
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viscous flux F given by 

F = idiwu-P-,{p)+P^{p^). (1.3) 

By adding and subtracting terms, we can rewrite the momentum equation as in in terms of F 
and Lo: 

pu= = F,^ + fiu;i'^ + pP - P^{p,),^ . (1.4) 
The decomposition ()1.4p also implies that 

/\F = dw{pu + pf -VP^{p^)). (1.5) 
We refer to Hoff [1] for a more detailed discussion of F. 

We now give a precise formulation of our results. First concerning the pressure P^ we assume 

that 

(1.6) there is X > such that P^{p) ^ Lpf for 7 > 1. 

We also fix a positive reference density poo and choose bounding densities < p < p satisfying 

P< Poo < P (1-7) 

and we define a positive number S such that 

5 = imin{|poo - p|, IPoo - (1-8) 
Concerning the external force /, we assume that 

(1.9) / = -pVV', where V £ H^iR^); 

(1.10) there is a constant C > such that \tp{x) \ + |V-0(a;)| < C and \Dl^{x)\ < |VV'(a;)| for x S M^; 

(1.11) lim i(;{x) = 0. 

|a;|— ^00 

Concerning the diffusion coefficients fi and ^, we assume that 

0<C< + (1-12) 

It follows that 

ip + (£.- N 

for p = 6 and consequently for some p > 6, which we now fix. 

We define p-y as mentioned at the beginning of this section. Given a positive constant densty poo , 
we say {pj,0) is a steady state solution to (jl.ip if G C^(R'^) and the following holds 

r VP^(p(a;)) - -p^(a;)VV'(x), 

<^ lim p(:E)=Peo. (1-14) 
By direct computation, pj has the following explicit form 

Poo cxp [--^^(a;)] if 7 = 1, 



(1.15) p^{x) = 



if 7 > 1. 



Concerning the initial data (po, uq) we assume that there is a positive number N, which may be 
arbitrarily large such that 

||uo||LP<iV, (1.16) 
and there is positive number 5 with d < 6 such that 

p + d < essinf po < ess sup po < p — d, (1-17) 
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We also write 

Cq^Wp^- p4l- + \\u4l.. (1.18) 

where is defined as in (1.15). 

Weak solutions are defined in the usual way; we say that (p, u) is a weak solution of (|l.ip - (ll.2p 
provided that (p - p^, pu) e C([0, oo); H^''- {R^}) with {p,u)\t=o = ipo,uo), Vtt £ L'^{M.^ x (0, oo)) for 
i > 0, and the following identities hold for times t2 > ti > and test functions ip having uniformly 
bounded support in x for t G [ti, t2]: 

t2 



p{x, ■)(p{x, ■)dx 



{pipt + pu ■ 'Vip)dxdt, (1-19) 



and 

t2 



{pw'){x, ■)(p{x, ■)dx 



tl 



[pw' ipt + pu-'u ■ Vip + P{p)ipxj]dxdt 

[pVu^ ■Vip+ {p- £.){div u)ipx.]dxdt (1.20) 
pifxjipdxdt. 



We use the usual notation for Holder seminorms: for u : M'^ M'" and a G (0, 1], 

/ \a \V{X2) - v{xi)\ 

(v) = sup — j j ; 

xi,x2e«3 \X2-Xi\°' 

and for t; : g C R3 X [0, oo) R™ and ai, ^2 e (0, 1], 



, vQi.as ^ \V{x2,t2) -v{Xi,ti)\ 

(^i,t,u^.t,>eQ \X2 - + \t2 - iih^ 



We denote the material derivative of a given function v hy v — Vt + Vw • u, and if X is a Banach 
space we will abbreviate by X when convenient. If J C [0, oo) is an interval, C^{I; X) will be the 
elements v G C(/; X) such that the distribution derivative vt G V{M.^ x int/) is realized as an element 
of C{r,X). Finally, if J7 C M'^ is a Lebesgue measurable subset in R'^, will be the corresponding 
volume of fl. 

The following is the main result of this paper, which is valid for the case when 7 = 1: 

Theorem 1.1 Let 7 = 1 and assume that the system parameters in (jl.ip satisfy the conditions in (1.6) 
and fl.QV ljl.lSp . Let positive numbers C, L, pac, p, p, S be given satisfying (|1.7p - (|1.8p and 



L log j < -C < C < L log j (1.21) 

Then given positive numbers N and d < 6, there are positive constants a, C, 6 depending on the 
parameters and assumptions in (1.6) and (1.9)- (|1.13p . on C, L, pac, p, p, R, on N and on a positive 
lower bound for d, such that if an initial data (po, uq) is given satisfying (fTTTell - ifLTSl) with 

Co < a, (1.22) 

then there is a solution {p,u) to (|l.ip - (|1.2p in the sense of (|1.19p - (|1.20p on all ofM.^ x [0, oo). The 
solution satisfies the following: 

p- p^, pueC{[0,oo);H-\R^)), (1.23) 

Vw G i^(M^ X (0,oo)), (1.24) 

u{-,t) G H\R^), t > 0, (1.25) 

uji-,t),F{-,t)eH\R^),t>0, (1.26) 
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(")rx[.,oo)<C(r)Co^ r>0, (1.27) 
where C{t) may depend additionally on a positive lower bound for t, 

p < p{x,t) <p a.e. onM.^ x[G,oo), (1.28) 

and 

sup / [{p - p^f + |up + a\Vu\^ + a^F^ + \Vuj\^)] dx 

+ / / [\Vu\^+a{\uf + \Vuj\^)+a^\Vu\^]dxds <CC^ (1.29) 
Jo Jr^ 

where a{t) = min{l, t}. Moreover, (p, u) — > (p-y, 0) as t oo in the sense that, for all ri G (2, oo] and 
ri e (2,oo), 

lim \\p{-,t)-~p^{-)\\L'-i+M-,t)\\Lr2 =0. (1.30) 

Using similar method, we also obtain parallel results for 7 > 1, with an extra assumption on the 
support of the external potential force -0: 

Theorem 1.2 Let 7 > 1 and assume that the system parameters in (jl.ip satisfy the conditions in (1.6) 
and fl.QV lll.lSp . Let positive numbers C, L, pao, p, p, S be given satisfying ()1.7p - (|1.8p and 

[pZo-' - -P'-'] <-c<c<^^ [pL-^ - P'~'] (1.31) 

Then given positive numbers N and d < 6, there are positive constants a, C, 6 depending on the 
parameters and assumptions in (1.6) and (1.9)- (|1.13p . on 7, C, L, poo, P, on N and on a positive 
lower bound for d, such that if an initial data (pQ,uo) is given satisfying (|1.16p - (jl.l8p and the small- 
energy assumption p.22p . and if 

|supp(V')| <a, (1.32) 
then there is a solution {p,u) to ()1. 1^ - 1)1. 2^ in the sense of (I1.19p - (|1.20p on all ofW^ x [0,cxd). The 
solution satisfies ()1.23p - ()1.27p with bounds p.28p - (|1.29p . and {p,u) {p-y,0) as t 00 in the sense 
of P^O)) for all ri e (2, 00] and ra £ (2, 00). 

We point out that for the case when 7=1, system (|l.ll) becomes simpler than those for 7 > 1. 
For example, in deriving a priori bounds for the term / / a\u\^dxds, we first rewrite the momentum 
equation in (|l.ip as follows 

pii + P^V{P-^P-^ - 1) + pP-^VP^{p-^P^ - p:,^P'y) - mAu -{£,- M)V(divu) = 0. (1.33) 
When 7 = 1, the third of the above reads 

pP^^VP-,{p-^P-, - p-^Py) = pP-^VP^{L -L)=0, 

so if we multiply (jl.33l) with aii and integrate, / / a\u\'^dxds can be bounded in terms of the L^ 

JoJm 

norms of (p — p, u) and the norms of u (modulo higher order terms in u). 

On the other hand, when 7 > 1, without any cancellation, we then have to estimate the term 

t 

crupP^^V Pj{p^^ Pj — p^^Pj)dxds which makes the analysis much more intricate and forces 

extra assumptions imposed on the external force. We will explain more later in subsequent sections. 

This paper is organized as follows. We begin the proofs of Theorem 1.1 and 1.2 in section 2 
with a number of a priori bounds for local-in-time smooth solutions. Since many of these estimates 
are rather long and technical, we omit those which are identical to or nearly identical to arguments 
given elsewhere in the literature. In section 3 we derive the necessary bounds for density by applying 
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the estimates in Theorem 3.1 and 3.2 in a maximum principle argument along particle trajectories of 
the velocity, making important use of the monotonicity of as described in (1.8). The small-energy 
assumption (|1.22p then enables us to close these arguments to show in Theorem 3.1 and 3.2 that both 
the pointwise bounds for density and the a priori energy bounds of Theorem 3.1 do hold as long as 
the smooth solution exists. Finally in section 4 we prove Theorem 1.1 and 1.2 by constructing weak 
solutions as limits of smooth solutions corresponding to mollified initial data. 

We make use of the following standard facts (see Ziemer [15] Theorem 2.1.4, Remark 2.4.3, and 
Theorem 2.4.4, also Ladyzhenskaya j7] section 1.4, for example). First, given r € [2,6] there is a 
constant C(r) such that for w G i/^(R'^), 

lklU.-(M3) < C{r) niVu;|!i^;-f;/^'-) . (1.34) 

Next, for any r G (1, cxd) there is a constant C(r) such that for w G W^'^{M.^), 

||w||l~(r3) < C(r)||w||^i,r(K3) (1.35) 

and 

(^)g3 <C(r)||Vu;|U.(R3), (1.36) 

where a = 1 — 3/r. If F is the fundamental solution of the Laplace operator on M"^, then there is a 
constant C such that for any / e L'^{R^) n L4(M3)^ 

nr., * /||l-(R3) < C (||/|U2(r3) + ||/|U4(R3)) . (1.37) 

Finally, there is a constant M such that for any v e i/^(]R.'^), 

J^g^dx<M/ \Vv\'dx. (1.38) 

3 (1 + \x\y 



2. Energy Estimates 

In this section we derive a priori bounds for smooth, local-in-time solutions of (jl.ip - (|1.2p whose 
densities are strictly posititve and bounded. Specifically, we fix a smooth solution [p — p^,u) of (|1.1|) - 
(|1.2I) on R'^ X [0, T] for some time T > with initial data {po — pj,uo). These bounds will depend only 
on the quantities Cq, C, L, pr^n p, p, N, d and will be independent of the initial regularity and the time 
of existence. 

We define a functional A{t) for a given such solution by 

A(t) = sup / [a\Vu\^ +a^{\u\^ + \VLu\^)]dx (2.1) 

0<S<t Jr3 

+ // [cr(|Mp + iVwp) +cr3|VMp] da;ds, 

where a{t) = min{l,t}, and we obtain the following a priori bound for A{t) under the assumptions 
that the initial energy Co in (jl.181) is small and that the density remains bounded above and below 
away from zero when 7 = 1: 

Theorem 2.1 Let 7 = 1. Assume that the system parameters in (jl.ip satisfy the conditions in (1.6) and 
(1.9)- (|1.13|) . Let positive numbers C, L, poo, p, p,S be given satisfying (|1.7|) - (|1.8p and ()1.21|) . Then given 
positive numbers N and d < 5, there are positive constants a, Af , 6 depending on the parameters and 
assumptions in (1.6) and (1.9)- (|1.13p . on C, L, pao, p, p, on N and on a positive lower bound for d, such 
that if {p — p-y, u) is a solution of ()l.ip - (|1.2p on X [0,T] with initial data {po — p^,uq) G H^{S.^)) 
satisfying (jl.l6l) - (|1.18p with Cq < a, and if 

P < P{x,t) < p on R3 X [0,T], 
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then 

A{T) < MC^. 



With an extra assumption on supp('0), we can also obtain the following estimate for 7 > 1: 

Theorem 2.2 Let 7 > 1. Assume that the system parameters in (II. 1|) satisfy the conditions in (1.6) 
and (1.9)- (|1.13|) . Let positive numbers C, L, p^o, p, p,S be given satisfying ()1.7p - ()1.8p and (|1.3ip . Then 
given positive numbers N and d < S, there are positive constants a,M,9 depending on the parameters 
and assumptions in (1.6) and (1.9V (I1.13I) . on 7, C, L, poo, P, on N and on a positive lower bound for 
d such that if |supp(?/')| < a and if {p — Pj^u) is a solution of (|l.ip - ()1.2p on x [0,T] with initial 
data {po - Pj,uo) € i/^(IR^)) satisfying (I1.16l) - (|1.18p with Cq < a and 

P < Pix,t) <p onM.^ X [0,r], 

then 

A{T) < MC^. 



The proof will be given in a sequence of lemmas in which we estimate a number of auxiliary 
functional. To describe these we first recall the definition (jl.lSp of p, which is an open condition, and 
which therefore allows us to choose q G (6,p) which also satisfies (|1.13p . Then for a given {p ~ p^,u), 
we define 



Hit) 



a^/^\\7u\^ + a^lVuHdxds 



D{t)= f f \p- p^\'^\SIP{p)\^dxds. (2.3) 

Also, we can readily see that assumptions (|1.2ip and (|1.3ip imply p^ is well-defined and 

P<Pi<P, (2.4) 

which will be crucial to the later analysis. 

For simplicity, we write a = a[t) = min{l,t}, P-y = P-y{p) and P-y = P-y(p-y) without further 
referring. 

We begin with the following L'^ energy estimate, which is valid for 7 > 1: 

Lemma 2.3 Assume that the system parameters in (|l.ip satisfy the conditions in (1.6) and fl.QV lll.lSp . 
Then if [p — p^, u) is a solution of (|l.ip - (|1.2p on x [0, T] with initial data (po — p^, uq) G H'^(R.^)) 
satisfying (|1.16p - (|1.18p . and if p, p G [p^p], then for 7 > 1, 

sup / {\p ~ p-yl'^ + p\u\'^)dx + [ [ \Vu\^dxdt < MCo- (2.5) 

0<t<T JR3 Jo 

Proof. We use V'0 — —p^^VP-y on the momentum equation to get 

pu + p{p^^\7P^ - p:,^VP^) - pAu - (C - m) V(div u) = 0. (2.6) 
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Multiply the above by u and integrate to obtain that for < t < T, 



\p\u\^dx 



dx 



-^0 



+ // [^i\Vu\^ + p){divuy\dxds ^ 0, 



(2.7) 



where the divergence of a matrix is taken row-wise. Next we define 

G{p) = n s-ip;(s)dsdr, 
then using the mass equation, the second term on the left side of (12.71) can be written as follows 

pu{p^^\7P^ - p:^^VP^)dxds 



1 



f [ pu-V \ f r-'^Pl^{r)dr] dxds 

oJR3 yjp^ J 



t r / rp \ ft /. f t 

pt I / r^^P!^{r)dr dxds 



'oJm.^ \-j P^ 
Putting the above into (|2.7p . the result follows 



G{p)sdxds — I G{p)dx 



□ 



Next we derive preliminary L'^ bounds for Vu and ii: 



Lemma 2.4 Assume that the hypotheses and notations o/ Lemma 2.3 are in force. Then for t G (0,T] 
and 7 = 1, 



sup a I |Vu|^dx+ / / ap\u\^dxds<M[Co+H], 
o<s<t Jr3 Jo- 



(2.8) 



while for 7 > 1, 



and 



sup cr^ / \u\'^dx+ / / a^\Vu\^dxds < M[Co + H] 

0<s<t JR3 J0JR3 



sup cr [ |Vupdx+ [ I ap\u\^dxds < M[Ca+ H + D] 
0<s<t Jr3 JoJm.3 



sup a I \u 

0<s<t 



a\^dx+ [ [ cr^\Vu\^dxds < M[Co + H]. 



(2.9) 



(2.10) 



(2.11) 



Proof. The proofs are nearly the same as those of (2.9) and (2.12) in |4|, except that is not 
necessarily a constant. We prove (j2.8p and (|2.10l) as an example. First, (|2.6I) can be rewritten as 
follows 

pii + P^V{P-,P-^ - 1) + pP-^WP-,{p-^P^ - p-^P^) - /iAu - (e - M)V(divu) = 0. (2.12) 
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We multiply (|2.12p by au and integrate to obtain 

/ / ap\u\^dxds — n / / aiidxds ~ ^ fi) / / auV{divu)dxds 
JoJr3 JoJR3 JoJr'' 
t . 

(juP^V{P^P:;^ - l)dxds 



aupP-^VP.y{p~^P^ - p~^P^)dxds = 0. 



(2.13) 



The second and the third term on the left side of (|2.13p can be bounded above by 

— T-cr I \SIu\'^dx — ^ ^ a I {divuf'dx + H. The forth term on the left side of p.lSp can be esti- 

mated as follows 



auP-y\7{P-yP:^^ - l)dxds 



< M 



a I \\/u\\p — p^\dx + I / IVuWp — p^\dxds 

t 

iVul' 



where the right side is bounded by M [Cq + A] by Lemma 2.3. Finally for the fifth term, we notice 
that for 7 = 1, 

p-^P., - p-^Py = L - L = 0, 

so that the term vanishes for the case when 7 = 1. On the other hand, for 7 > 1, using the definition 
((23| of L), we have 



aupP:^^V P^{p^^ P^ - p~^P^)dxds 







[fl 


a\u\'^dxds 




3 



<M / a\u\(j){-f)\p ~ p-f\\Vp^\dxds < 
Therefore jMl) and ((2T0l) follow. 

The following auxiliary estimates will be applied to bound the the functional H: 

Lemma 2.5 Assume that the hypotheses and notations 0/ Lemma 2.3 are in force. Then for 
< i < 1 AT and 7 > 1, 

sup / \u\idx+ ff [\u\i-^\Vu\^ + \u\''-^\Vi\u\^)\^]dxds 

<S<t JR3 J0JR3 

< M 



□ 



0< 



(2.14) 



P-1 q-2 

C^-'N^ +Co 



Proof. We multiply (|2.12p by ^ and integrate to obtain that 



\u\'^dx 



< M 



{P^ ~ P.y)div{\u\''''^u)dxds 



J0JR3 



(2.15) 



'0 



\'i-^\S/p^\\P-y -P^\dxds 



f f pP-^ip-^P^ - p-^P^)\u\'i-^udxds + f \uo\'^dx 
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The last term on the right side of (|2.15p is bounded by ^ N p-^ , and the first and second term can 

be bounded in terms of sup / \p — pj\'^dx. For the remaining term we have 

o<s<t Jm 



pp-\p-^P^^p-^P^)\u\'i-^udxds 

^ [ [ \^P7\Hl)\P - Pf\H'^~^dxds 
JqJr^ 



< M 



sup 

0<s<t , 



1 





Q 


sup / 


lul^dx 


/R3 




.0<s<t Ji 


.3 



Using the above on (|2.15p and absorbing terms, the result follows. 



□ 



We obtain bounds for u.uj and F in W^'"^, these being required for the derivation of estimates 
for the auxiliary functionals H and D: 

Lemma 2.6 Let 7 > 1 and assume that the hypotheses and notations 0/ Lemma 2.3 are in foree.Then 
for r G (1, 00) and t G (0, T], 



\Nu{;t)\\Lr < M [m;t)\\Lr + M;t)\\Lr + \\{P, - P-,){;t)\\Lr 

The constant M in (|2.16p may depend additionally on r. 



(2.16) 



Proof. We have from the definition ()1.3p of F that 

^Au^ =F,^+^uji^^ + {P^-P,),^. 
Differentiating and taking the Fourier transform we then obtain 

^uliy,t) = y2yiFiy,t)+^f^:^^iy^t) + yj^iP~P^)iy,t). 



(2.17) 



The result p.l6p then follows immediately from the Marcinkiewicz multiplier theorem (Stein [12] , pg- 
96). □ 

Lemma 2.7 Assume that the hypotheses and notations 0/ Lemma 2.3 are in force. Then for 7 > 1, 



\\VF(;t)\\L- + \W^{-ML- 



< M 



\\u{;t)\\L2 + \\Vu{-,t)\\L2 + \\\\7P\\p-p^{;t)\\L2 



while for 7 = 1, 



IV ( |j {■,t)\\L2 + < M[\\ui;t)\\L2 + \\Vu{;t)\\L2] . 



Proof For ([STS]) . using (fTS]) we have 

AF = ipu^),^ + [iP^),^p-\p^ - P)h,. 



(2.18) 



(2.19) 



(2.20) 



Similar to p.l7p . we differentiate and take Fourier transform on ()2.20p and apply the multiplier 
theorem, we get the bounds for VF and similarly for divw. To show ()2.19p . using (|1.3p and the 
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definition (|1.4p of i^, for 7 = 1, 



P 1 ^ 



Lp: 

P ( iO 



-^div u — pp ^ + 1 
Lp 



P ( ;^-^P ) + Lp^^ipp ^ - 1) 



and so, 



PJ X, L\P 



|V(|)(-,t)|U. <Af 



H ^divw H ^ — uj, 

P P 



\\u{;t)\\L2+\\Vui;t)\\L2 + \\V [j]{;t)\\L2 



Using similar method, we obtain 



IV ( ^ ) i;t)\\L2 < M[\\u{;t)\\L2 + \\Vu{-,t)\\L2]. 



Therefore ((2l9| follows from ([2:221) and ^22il- 



(2.21) 
(2.22) 

(2.23) 
□ 



Next we derive a bound for the functional D defined above in (12.3 



Lemma 2.8 Assume that the hypotheses and notations 0/ Lemma 2.3 are in force. Then for "f — 1 and 

te{0,T], 



D{t) < M[Co + A{t)]. 
For "f > 1, (j2.24p also holds when |supp(?/')| is sufficiently small. 



(2.24) 



Proof. In view of Lemma 2.3, it suffices to consider t > 1. We first consider the case when 7=1. 
Using the mass equation. 



D ~ F 



(2.25) 



so that we multiply the above by (p — p^)sgn(p — p-y)|VPp and integrate to obtain 

Jr3 JiJr^ 
< M 



.F . 



< M 



+ \u\^){l + \x\)~''dxds 
iJk? P 

VP. 



|Vwr)dxds 



(2.26) 



where the last line follows by ([QS]) . Using ([23]) and ([2TT9|) . the right side of (|2?26l) is bounded by 
M [Co + A], and hence (PTMl) holds for 7 = 1. 
For 7 > 1, using (fOS]) . ([21^ and 



p- p^Y\yP\^dxds < M 



< M 



< M 



Co 

Co- 

Co 



|Fp + \u\'^)\VP\'^dxds 



Vu\^dxds+ I I \F\^\VP\^dxds 

1JR'> Jl 
t , 

\F?\\/P\'^dxds 

1 •/ SUpp('0) 



(2.27) 
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where the last inequality follows from the fact that |VP| = /5|VV'|. Assume that |supp('(/')| < oo, the 
last term on the right side of (|2.27p can be estimated as follows 



1 J supp(^) 



|Fp|VPpdxds < Isupp(V') 

< |supp(V')|^ 

< |supp(V')|^ 

< |supp(V')|^ 



\F\^\VPVdx 



/supp(?/;) 

\F?\VP?dx 



ds 



\F\'^\DlP? 



ds 



\VuY + \ypY\p~ pY)dxds 



\yPY\p - pY)dxds + M(Co + A) 



\F\^\VP\^dxds + M 

\F\^\VP\^dxds + M 
hence if |supp('!/')| is sufhciently small, we have 

// \p- p^\'^\^P\'^dxds <M[Cq+ A] 

and ([2:241) follows. 

We can now obtain the required estimates for the functionals H defined above in (12.2 



□ 



Lemma 2.9 Assume that the hypotheses and notations o/ Lemma 2.3 are in force. Then for "f — 1 and 

\21 



H{t) < M [Co + Cl + A{tf 
For "f > 1, (j2.28l) also holds when |supp(?/')| is sufficiently small. 



(2.28) 



Proof. We first consider the case for 7 = 1. In view of the definition (|2.2p of H, we bound only the 
o-^/^IVmI^ + cr^|Vu|'*] dxds. 



term 



The remaining term 



dxds 



First we have from Lemma 2.6 that 



a^\Vu\'^dxds < M 



is bounded exactly as in HofF [4] pp. 29-32. 



(73(|p-p^|4+|F|4+|w|4) dxds 



(2.29) 



Applying ()1.34p we can bound the second term on the right by 



F f f F 

(T^\ — \'^dxds<[ sup / a\ — \^dx 

0Jr3 P \0<s<tjR3 P 



Now from the definition of F and Lemma 2.3, 



IV ( 1^ ] \^dx 



crlV 



dxds 







i sup / 


" a\-\''dx] 


VO<s<i Jl 


13 P J 



Also, from ^J^, 



o-|V ( ^ 1 [■'dxds < M 

0JR3 \ P 



cr + |VmP) dxds 



<M[A + Co] 
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and 



Thus 



sup / + |Vup) 

0<s<tjR3 



< MA. 



Applying the resuhs of Lemma 2.7 in a similar way, we obtain that 



(2.30) 



p' I I cr^l^l^dsds < ( sup / a\'^Ydx\ ( |V ( ^ ) l^dx 



OJ 



< M[Co + A]' 



(2.31) 



For the first term on the right of p.29p . multiply (|2.32p by (T^\p — p-y|'^sgn(p — p^), integrate and use 



a^\p — p-yl'^dxds < M 



< M 



Co + / / a^QF"^ + \u\'^)dxds 

2 



Co + (Co + AY 



(2.32) 



Hence we obtain a bound for j j (7'^\'Vu\'^dxds from ()2.30p . p.3ip and (|2.32p . Bounds for the term 
/ / cr'^^^l'Vul^dxds are obtained in a similar way, thereby proving ()2.28p for 7 = 1. 



For 7 > 1, using ((2TT8)) . 

a^\F\^dxds<[ sup / a\F\'^dx 

0JR3 \0<s<t 



< (Co+A)^ 



alVFrdx ds 



\JR3 
2 , |V7„.|2 



\JR3 

1 , , 3 



a{\uY + \VuY + \p- pY\VPYdx] ds 



< MiCo + A)2 (A's + A^Co + D) 
<M{Co+A + Df, 



and similarly. 



a^\p- p^\'^dxds+ I I a^\uj\^dxds<M{CQ + A + Df. 



So when |supp(V')| is sufficiently small, we can apply p.24p to conclude that 



0"'R=' 



a'^lVuYdxds < M[Co + (Co + Af 



which proves (|2.28p for 7 > 1. 



□ 



Proof of Theorem 2.1 and 2.2: Theorem 2.1-2.2 now follows immediately from the bounds (|2.5p . (|2.8p . 
((2^ . (I2.14p . (|2.24l) and (I2.28p . and the fact that the functional A is continuous in time. □ 
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3. Pointwise bounds for the density 

In this section we derive pointwise bounds for the density p for 7 > 1, bounds which are independent 
both of time and of initial smoothness. This will then close the estimates of Theorem 2.1-2.2 to give 
an uncontingent estimate for the functional A defined in (j2.1[) . The result is as follows: 

Theorem 3.1 Let "f — 1. Assume that the system parameters in (|1.1|) satisfy the conditions in (1.6) and 
(1.9)- (I1.13I) . Let positive numbers C, L, poa, p, p,d be given satisfying (jl.7l) - (|1.8p and (|1.21l) . Then given 
positive numbers N and d < S, there are positive constants a, M, 6 depending on the parameters and 
assumptions in (1.6) and (1.9V (|1.13p . on C, L, pao, p, p, on N and on a positive lower bound for d, such 
that if {p — p-y, u) is a solution of (|l.ip - (|1.2p on X [0,T] with initial data (po ~ P-y^uo) G H (M. )) 
satisfying (jl.l6p - (|1.18l) with Cq < a, then in fact 

P < P{x,t) < p on R3 X [0,r], 

and 

A{T) < MC^. 

Theorem 3.2 Let 7 > 1. Assume that the system parameters in (|l.ip satisfy the conditions in (1.6) 
and (1.9)- (jl.l3p . Let positive numbers C, L, poa^ p^ p,S be given satisfying (ll.7p - (ll.8p and (|1.31l) . Then 
given positive numbers N and d < 5, there are positive constants a, M, 9 depending on the parameters 
and assumptions in (1.6) and (1.9)- (|1.13|) . on 7, C, L, p^o, P, P, on N and on a positive lower bound for 
d such that if |supp(V')| < a and if (p — p-y^u) is a solution of p.ip - (|1.2|) on M.^ x [0,r] with initial 
data (po — PjtUq) € i/'^(IR'^)) satisfying ()1.16|) - (|1.18p with Cq < a, then in fact 

P < p{x,t) <p onR^ X [0,r], 

and 

A{T) < MC^. 

Proof. We prove Theorem 3.1 as an example, the proof of Theorem 3.2 is just similar. Most of the 
details are reminiscent of those in Hoff-Suen ,13j section 3. First, we choose positive numbers b and b' 
satisfying 

p<b<p + d<p-d<b'<p. 

Recall that po takes values in [p + d, p — d], so that p g [p,p] on R'^ x [0,t] for some positive t. It 
then follows from Theorem 2.1 that A{t) < MC^, where M is now fixed. We shall show that if Cq 
is further restricted, then in fact 6 < p < 5' on R"^ x [0, r], and so by a simple open-closed argument 
that 6 < p < 6' on all of R^ x [0, T], we have A(T) < MCq as well. We shall prove the required upper 
bound, the proof of the lower bound being similar. 

Fix y gR^ and define the corresponding particle path x{t) by 

f x{t) =u{x{t),t) 
{ x{0) = y. 

Suppose that there is a time ti < t such that p{x{ti),ti) = b'. We may take ti minimal and then 
choose ^0 < ti maximal such that p(a;(io), to) = p — d. Thus p{x{t), t) e [(p — d), b'] for t G [to, ^i]- We 
consider two cases: 

Case V.to<ti<TM. 

We have from the definition (II. 3p of F and the mass equation that 

p^\\og p{x{t),t)] + P{p{x{t),t)) ~ P^{p^{x{t)) = -F{x{t),t). 
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Integrating from to to ti and and abbreviating p(x(t),t) by p{t), etc., we then obtain 



/i[log(&') - \og{p -d)] + I ' [P^{s) - P^{s)]ds = - r F{s)ds. 



We shall show that 



F{s)ds < MC^ 



(3.1) 



(3.2) 



for a constant M which depends on the same quantities as the M from Theorem 2.1 (which has been 
fixed). If so, then from p.l|) . 



M[log(6') - log(p - d)] < A^Co^ 



(3.3) 



where the last inequality holds because P^(s) — P-y(s) is nonnegative on [to, ti]. But p. 31) cannot hold if 
Co is small depending on M, and p — d. Stipulating this smallness condition, we therefore conclude 
that there is no time ti such that p{ti) — p{x{ti),ti) — b' . Since y e M"^ was arbitrary, it follows that 
p < b' on R'^ X [0, r], as claimed. The proof that 6 < p is similar. 

To prove (|3.2I) we let F be the fundamental solution of the Laplace operator in R'^ and apply 
(11.51) to write 



F{s)ds = 



to 



^x,{x{s) - y)pu^{y,s)dyds 

^x.ixis) ~ y) (P-y):r,p7^(P7 - P) 



dyds. 



The first integral on the right here is bounded exactly as in Lemma 4.2 of HofF [4]: 

{x{t) - y)pu^ {y, t)dyds 

< \\T.,^*{pu^){-M)\\L^ + \\^x,*{pu'){-,to)\\L^ 

f f T^^^,{x{s) ~y)[u\{x{s),s) -u''{y,s)]{pu^){y,s)dydi 
to Jr^ 

< MCl 

and the second integral on the right side of p.4p can be bounded as follows 



^x,{x{s) -y)[{P^)^^p^^{p^ - p)]dyd. 



<M \y\-^\p-p^\^dyds 





4 




/ \y\ '^dy 




/ 1^ 


l\v\<^ 




-'|y|<i 



\p~Pi\ dy 



< M 



+ M 
< MCl 

where the last inequality follows from Theorem 2.1. Thus p.2p is proved. 
Case 2: 1 < to < ti. 



ds 







2 




2 


f 


[ \y\-'dy 




/ \P-Pf\'^dy 


ds 


Jo 


y\v\>i 




Av\>i 





(3.4) 
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Again by the mass equation and the definition (|1.3p of F, 



-{p{t) - p^{t)) + p-^ p{t){P^{t) - P^(t)) 



= -p-^p{t)F{t) ~ p-^u(t) -Vp^it). 
:{p{t) - p,{t)r + p-'g{i)p{t){p{t) - p,{i)f 



Muhiplying by {p{t) - p^{t)f we get 

= -p-^p{t){p{t) - p^{t)fF{t) ~ p-^u{t) ■ Vp^{t){p{t) - p^{t))\ 
where g{t) = (P^(i) - P^{t)){p{t) - p^{t))-'^ > on [io,ii]- We integrate to obtain 
{p{h) ^ p^{h)f - {p{to) - p^,{h)f 



< M 



\F{;s)\\l + \\u{;t)\\to] ds. 



to 



We shall show that 



so that from (13.61 



M 



\F{;s)\\l + \\u{;t)\\U ds.<MCl 



to 



(3.5) 



(3.6) 



(3.7) 



Q<{h' ~pf ~{p~d~ pf < MCl 

This cannot hold if Co is sufficiently small, however, so that, as in Case 1, there is no time ti such that 
p{ti) = p{x{ti),ti) = 6'. Since y G M'^ was arbitrary, it follows that p < 6' on x [0,r], as claimed. 

To prove p.7p . it suffices to consider the term M f \\F {■ , s)\\'^ds . We apply (|1.37p to get 



to 



tl r 



to 



\\F{-,s)\\l,ds< / \\pu{-,s)\\l, + \\VP^~p-\~p^-p){;S)\\l, 

r*i r 

I to 



\\pu{;s)\\l. + \\VP^p-\p^- p){;s) 



ds 



ds, 



(3.8) 



and the right side of the above is readily seen to be bounded by MCq . Thus (|3.7p is proved. □ 



4. Global Existence of Weak Solutions: Proof of Theorem 1.1-1.2 

In this section we complete the proof of Theorem 1.1 and 1.2 by constructing weak solutions as limits 
of smooth solutions. Specifically we fix the constants a and M defined in Theorems 2.1-2.2 and 3.1-3.2, 
we let initial data {po,uo) be given satisfying the hypotheses (|1.16p - (|1.18p and (|1.22p of Theorem 1.1 
and Theorem 1.2, and we take (pq, Uq) to be smooth approximate initial data obtained by convolving 
(POjUq) with a standard mollifying kernel of width ?/ > 0. We apply the local existence results (see 
Nash [H] or Tani [M]) to show that there is a smooth local solution {p^ , u^) of (|l.ip - (|1.2p with initial 
data (pq,Mq), defined up to a positive time T, which may depend on 7]. The a priori estimates of 
Theorem 3.1-3.2 then apply to show that 

A{t) < MC^ and p < p'^{x, t) < p, (4.1) 

where A{t) is defined by p.ip but with {p,u) replaced by {p^,u^). By standard arguments together 
with the bounds (|4.ip . [p^^u^) exists and satisfies (|l.ip - (|1.2l) for all time. 

Those bounds in (14. ip will provide the compactness needed to extract the desired solution [p, u) 
in the limit as 77 — > 0. We begin by proving uniform Holder continuity of the families {u^} away from 
t = 0: 
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Lemma 4.1 Given t > there is a constant C — C{t) such that, for all rj > 0, 

("''(■-0)M[.,oo)<^^MCo- (4.2) 
Proof. The proof is exactly as in HofF [4j, pg. 33 and pg. 41-42. □ 

Compactness of the approximate solutions {p^,u'^) now follows: 

Lemma 4.2 There is a sequence 77^ — > and functions u and p such that as k ^ 00, 

u^'' — > u uniformly on compact sets in M.^ x (0, cx)); (4-3) 

Vu'"= (■, i), Vw"^ {■,t)^ Vu(-, t), Vuj{-,t) (4.4) 

weakly in L^(R^) for all t > 0; 

a^u^^ ^a^VvP'' ^ a^iiaiVu (4.5) 

weakly in i^(M^ x [0,oo)); and 

pn-{.,t)^pi;t) (4.6) 

strongly in Lf^^{M!^) for every t > 0. 



Proof. The uniform convergence (|4.3p follows from Lemma 4.1 via a diagonal process. The statements 
in (14.41) and (|4.5p then follow for this same sequence from (|4.1I) and elementary considerations based 
on the equality of weak-i^ derivatives and distribution derivatives. The convergence of approximate 
densities (|4.6p for a further subsequence requires a more involved argument, given in Lions [8 pp. 
21-23 and extended by Feireisl [2], pp. 63-64 and 118-127. □ 



Proof of Theorem 1.1-1.2: We only prove Theorem 1.1 since Theorem 1.2 can be proved in a similar 
way. It is clear that the limiting functions {p,u) of Lemma 4.2 inherit the bounds in ()1.24|) - ()1.29|) 
from (|4.ip and (|4.2I) (but notice that no statements are made in (|1.28l) concerning u{-,t)). It is also 
clear from the modes of convergence described in Lemma 4.2 that {p,u) satisfies the weak forms 
(|1.19l) - (|1.20l) of the differential equations in (|l.ip as well as the initial condition (jl.2p . The continuity 
statement ()1.23|) then follows easily from these weak forms together with the bounds in (|1.29p . 

It remains to show (|1.30p . By the same argument as in Hoff [4] pp. 44-47, we have for r2 G (2, 00), 

lim \\p{.,t)- p^{-)\\L^+\\ui;t)\\L^, =0. 

i— >oo 

So using (fL29l) . for ri e (2, 00), 

hm \\p{.,t) - p^{■)\\L'■^ <CC'\im ||p(.,i)-p^(-)||L~ =0. 

i— foo i— >-oo 

This completes the proof of Theorem 1.1. □ 
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